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(N ! Abstract 

A general critical condition of subharmonic oscillation in terms of the loop gain is derived. Many closed-form 
critical conditions for various control schemes in terms of converter parameters are also derived. Some previously 
known critical conditions become special cases in the generalized framework. Given an arbitrary control scheme, a 
systematic procedure is proposed to derive the critical condition for that control scheme. Different control schemes 
share similar forms of critical conditions. For example, both V 2 control and voltage mode control have the same 
vq [ form of critical condition. A peculiar phenomenon in average current mode control where subharmonic oscillation 

{SJ ■ occurs in a window value of pole can be explained by the derived critical condition. A ripple amplitude index to 

predict subharmonic oscillation proposed in the past research has limited application and is shown invalid for a 
converter with a large pole. 
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I. Introduction 

J^i | Average continuous -time models are generally applied to analyze DC-DC converters |QX The average model can 
t-h ' predict some types of instabilities, but it generally cannot predict subharmonic oscillation [|2| which is associated 
^O , with a discrete-time pole at -1. By considering the sampling effects and increasing the system dimension, improved 
; | average models can predict the occurrence of subharmonic oscillation in some cases 01, EJ- However, even with 
£D ' positive gain or phase margins in the improved average model, the subharmonic oscillation may still occur 0. 
^vq This paper shows a systematic approach to accurately derive the critical condition. 

*— I \ Some subharmonic oscillation critical conditions for particular control schemes have been known. These control 
J> ■ schemes are current mode control (CMC) |[T] (where the subharmonic oscillation occurs at duty cycle D = 1/2), 
voltage mode control (VMC) without considering the effect of equivalent series resistance (ESR) J6], and V 2 control 
0. These conditions are generally obtained case by case, not systematically. 
This paper tries to answer the following questions: 

1) Is there a general closed-form critical condition that directly leads to these known conditions, and these 
known conditions become special cases in the generalized framework? 

2) There exists a peculiar phenomenon in average current mode control such that subharmonic oscillation occurs 
in a window value of pole [8]. Is there a general theory to explain this phenomenon? Does it also exist in 
other cases? 

3) It is hypothesized in @] that the signal ripple amplitude can predict the occurrence of subharmonic oscillation. 
Does it have limited application only in some particular cases? 

4) Some control schemes (such as V 2 control, VMC and CMC) may seem different, but do they share the same 
or similar form of critical condition? 

5) Given an arbitrary control scheme, is there a systematic method to derive the critical condition for that control 
scheme? 

The answers to all of these questions will be shown to be affirmative. 
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Figure 1. An illustrative critical boundary in the parameter space (v s ,R). 




Figure 2. A VMC buck converter with a compensator G c (s). 



Like the describing function approach J9], harmonic balance ifTOl . ifTTI is a tool to analyze a nonlinear system. 
Based on harmonic balance, a critical condition of subharmonic oscillation in the buck converter is obtained in 
CGI, lfl3l . lfT4l . Since most converter designers are familiar with the loop gain analysis, here the critical condition 
in H3, fl3l . 11141 is expressed in terms of the loop gain. Based on this critical condition, many closed-form critical 
conditions for various control schemes are obtained. To determine the loop gain for CMC, there exists different 
views about the PWM modulator gain lTl3Tl . For stability analysis, it will be shown that the PWM modulator gains 
for CMC and VMC are the same. 

Note that, here, the critical condition is expressed in the converter parameter space. For example, given a source 
voltage v s and load resistance R, an illustrative critical boundary in the parameter space (v s , R) is shown in Fig. [Q 
The critical conditions defines the subharmonic oscillation boundary in the parameter space to separate stable and 
unstable regions. When a converter parameter crosses the critical value, the stability (or instability) changes. Critical 
conditions in closed-forms greatly facilitate the converter design, because the quantitative effect of each relevant 
converter parameter can be clearly seen. 

The remainder of the paper is organized as follows. In Sections JI] and [Till, the operation of the buck converter 
and harmonic balance analysis are briefly reviewed. In Section ITVl a general critical condition in terms of loop gain 
is derived. In Section |VJ the critical conditions for typical loop gains are derived. In Sections IVlllIXI the proposed 
approach is systematically applied to various control schemes. Conclusions are collected in Section |X] 



II. Brief Review of Buck Converter Operation 

Consider a VMC buck converter shown in Fig. [2] where v s is the source voltage, Vd is the voltage across the diode, 
v is the output voltage, v r is the reference voltage, y is the compensator output signal, and G c (s) = —y(s)/v (s) 
is the compensator transfer function. Denote the ramp slope as m a = h(d). Denote the cycle period as T, the 
switching frequency as f s = l/T and let u s := 2irf s . 

Similarly, a CMC buck converter is shown in Fig. [3] where a control signal i c controls the (peak) inductor current 
il, and y = ic — ib (equivalently, G c (s) = 1). Note that in Fig. [3] the circuit is arranged in a way to have a similar 
output signal y as in VMC. 

The operation of the converter is as follows. Within a cycle period T, the dynamics is switched between two 
stages, S± and 52. Switching occurs when the compensator output y intersects with the ramp signal h(t) := 
Vi + {Vh — Vi)(y mod 1), where h(t) varies from a low value Vi to a high value T4- Denote the ramp amplitude 
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Figure 3. A CMC buck converter. 
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Figure 4. An equivalent nonlinear switching model. 



as V m = Vh — Vi. In the trailing-edge modulation (TEM) |0, Vd = f s in stage Si and i>d = in stage S^- The 
waveform of Vd{t) is a square wave with a duty cycle D. In the leading-edge modulation (LEM), Vd = in stage 
S\ and Urf = v s in stage S^- 77z/s paper focuses on TEM only. For LEM, see ^72l/ , ^731/ , 041/ . The controlled 
buck converter is equivalent to a nonlinear switching model shown in Fig. [4] Note that v r does not affect the loop 
stability, and it adjusts the DC offset of h(t) in Fig. ES Let the loop gain be T(s). From Fig. HI T(s) = v s G(s)/V m , 
and the model in Fig. |4] can be normalized as in Fig. |5] 

III. Brief Review of Harmonic Balance Analysis 

A brief review of harmonic balance analysis based on lILTI . lTT3l . lfl4l is presented. The square-wave Vd(t) can be 
represented by Fourier series (harmonics). The signal path in the control loop has two parts: from y to Vd through a 
nonlinear PWM modulator, and from Vd to y through a linear transfer function G(s) := G p (s)G c (s), where G p (s) 
is the power stage transfer function and its representation depends on what signal is fed to the compensator. 

In VMC, v is fed to the compensator. Let p = R/(R + R c ). For R c = 0, p = 1. From [fl] p. 470], the power 
stage Vd-to-v transfer function is 

r<( ,_ sRcC + 1 
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In CMC, %l is fed to the compensator (with G c (s) 
function is 
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Let x°(t) be the T-periodic solution of the converter. Let y°(t) be the corresponding T-periodic compensator 
output signal. The intersection of h(t) with y°(t) determines the duty cycle and hence the waveform of Vd{t). By 
"balancing" the equation y°(t) = h(t) (written in Fourier series form) at the switching instants, conditions for 
existence of periodic solutions and subharmonic oscillation can be derived. Let Re denote taking the real part of 
a complex number. Based on lfT2l . |[T3l , lfT4"l . the subharmonic oscillation occurs when 
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This critical condition is valid for a general switching system shown in Fig. |4] It will be shown later that this 
critical condition can be expressed in a closed-form related to the hyperbolic function csch . 
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Figure 5. An equivalent normalized nonlinear switching model. 

IV. Critical Condition in Terms of Loop Gain 

Since T(s) = v s G(s)/V m , the critical condition © directly leads to the following result. 

Consider a closed-loop buck converter with a loop gain T(s). For V m ^ 0, the critical condition of subharmonic 
oscillation is 
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Note that (@]) is an expression of convenience. For V m = 0, the loop gain v s G{s)/V m would be infinite. In that 
case, the equivalent critical condition © is used. Also note that (0]) is valid for both VMC and CMC by assuming 
the same PWM modulator gain to be 1/V m . Since both © and © are exact critical conditions, they can be used 
as benchmarks to determine the accuracy of other critical conditions. 

Generally, the power stage has an order of at least two (associated with L and C). If the compensator has an 
order of three, such as the type-Ill compensator [161, the loop gain T(s) has an order of five. However, there are 
ways to simplify the analysis without losing the accuracy. 

First, the compensator is generally designed to cancel some poles or zeros of the power stage, and the order of 
T(s) is reduced. Second, from ©, one sees that the shape of the Bode plot for the frequency smaller than uj s /2 
is irrelevant. One can use a simplified high-frequency form of T(s) for stability analysis. Third, by decomposing 
the loop gain into partial fractions, many new closed-form critical conditions in terms of the converter parameters 
can be obtained. 

Define an "F-transform" of T(s) as 



F[T(s)] := 2Re 
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Note that J r [T(s)] is a function of many converter parameters. Here, this function is called an L-plot, denoted 
as £. For example, let it be a function of D and it becomes £(D). Then, the critical condition © becomes 
C := T[T(s)] = 1. The critical condition itself does not tell which side of the critical boundary will be the stable 
region. Generally for a converter, the region with C < 1 is stable. The F-transforms of typical loop gain functions 
are presented next. 



V. "F-Transforms" of Typical Loop Gain Functions 

The loop gain is generally designed to have sufficient gain and phase margins. For a stable converter, the phase 
of T(juj) is less than 180° at the crossover frequency uj c . One can focus only those loop gains of first or second 
orders. 

The loop gain can be further decomposed into a combination of partial fractions. Only partial fractions of first 
orders are considered. Similar analysis can be applied to partial fractions of second orders. 

Let uj p and uo z be the pole and zero of T(s). Let p = ui p /oj s and z = uj z /uj s . The F-transform of the fraction 
l/(s + Up) will be the building block to derive the F-transforms of other loop gains. 



Table I 
F-TRANSFORM OF TYPICAL LOOP GAIN T(s). 



Case T(s) J-[T(s)] (note: p = uj p /oj s and z = uj z /uj s ) 

~ =*= ±a(D,p) = JR5) - Ql (Z3)p + C (Ap)) 
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:= a(Z>,p) 

where the proof of © can be obtained by looking up from a handbook of mathematical formulas or checked 

by a simple computer program, and csch is a hyperbolic function. Using Taylor series expansion, let a(D,p) = 

00 

£ {-l) k a k (D)p k . Using the L'Hospital's rule, one has a (D) = ir(2D - 1) and cti(D) = ir 2 (2D 2 - 2D + 1). 

fc=0 

A plot of a(D,p) is shown in Fig. [6] One sees that a(D,p) is close to zero for D > 0.5 and p > 1. The right 

straight-line edge in Fig. [6] is a(D,0) = ao(D) = ir(2D — 1). 

00 

Let the correction term be c(D,p) = Yl (—^-) k< ^k{D)p k . One has 

k=2 

c(D,p) = a(D,p)-a (D) + ai (D)p (7) 

A plot of c(D,p) is shown in Fig. [7J For a large p, since a(D,p) is small, c(D,p) increases almost linearly as a 
function of p shown in Fig. [7J From Fig. |7J the correction term c(D,p) is significant only if p > 0.1 (equivalently, 
Lo p > 0.1lo s ). For p < 0.1, c(D,p) can be ignored. 

The F-transforms of other loop gains are shown in Table Q] For each case C1-C9, the F-transform can be derived 
or proved (by simple algebra) in three ways. First, follow the definition of the F-transform as in ©• Second, 
by decomposing T(s) into a combination of fractions 1/s, l/(s + u> p ), etc, T[T(s)] is a combination of F[l/s], 
J r [l/(s + uj p )], etc. For example, 

• Use of C\ and C2 leads to C5 and C§; 

• Use of C2 and Cg leads to C7; 

• Use of Ci, C2 and Cg leads to Cg\ and 

• Use of C\ and the fact that F[l] = — 1 leads to C4. 

Third, each case is a special/general case of other cases. For example, by setting co p — > 0, C\ leads to C2', C§ leads 
to Cq; and C§ leads to C-j. By setting w 2 — > 00, C7 leads to Cq. One sees that Ci is a building block for other cases 
because they have similar terms as C\. 
Remarks: 

(a) All of the transforms in Table U are exact. No approximation is assumed. 

(b) There is no correction term c(D,p) for C2, Cq, or C-j. All other cases have a correction term c(D,p), which 
is small and can be ignored if p < 0.1 as discussed above. 

The critical conditions for various control schemes are readily derived next. 




Figure 6. Plot of a(D,p) for case &. 
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Figure 7. Plot of the correction term c(D,p). 

VI. CMC: CaseC 2 

From Fig. [3] T(s) = G p (s)v s /V m . From ©, at high frequency, G p (s) w 1/Ls. One has T(s) « t; s /LV^,,s, 
which is of case C2 in Table HI Its F-transform leads to the well-known critical condition for CMC ID, ||5l 

v s a (D) 
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Without the compensation ramp (m a = 0), the critical point is D = 1/2. 



VII. Proportional Voltage Mode Control (PVMC): Cases C7, C%, C$, and C3 

In proportional voltage mode control (PVMC), let the voltage loop have a proportional feedback gain k p . One 
has G c (s) = k p and y = k p (v r — v ). Five cases are considered to see the effects of R c and R on the subharmonic 
oscillation. 
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Figure 8. A buck converter under CF-PVR. 



A. PVMC with R c / 0: Case C 7 

From CO, at high frequency, G p (s) ss p{sR c C + l)/V m LCs 2 , and 
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which is of case C-j with to z = 1/R C C. From Table H the critical condition is 

v 8 kp P T 2 2RcC _ 2 _ 2D + 1)] = 1 

which is a weighted combination of 2D — 1 (a condition for CMC as seen in (©) and 2D 2 — 2D + 1. The term 
2D — 1 can be ignored only if R C C <C T. Therefore, in PVMC with a large R c , the critical condition has a term 
2D - 1 like CMC. 



B. V 2 Control: Case Cj, Same as PVMC 

The constant-frequency peak voltage regulator (CF-PVR, as shown in Fig. [8), a type of V 2 control Q, is proved 
below to be a special case of PVMC in terms of the critical condition. In CF-PVR, the output voltage is sensed 
(through a voltage divider with a gain of k p ), added with a stabilization ramp h(t), and compared with a reference 
signal v r to determine the duty cycle. One has G c (s) = k p , same as PVMC. Therefore, they have the same form 
of critical condition (fTTT ). 

For the buck converter, v = Dv s . To avoid subharmonic oscillation, (ITTb is rearranged as (inequality) 
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k p pv R c 

which agrees with [7. Eq. 7] and shows the required ramp slope m a to stabilize the subharmonic oscillation. 

Note that the critical condition (fT2l has two terms. The second term is related to the output capacitor C. For 
very large C, the second term can be ignored, and the critical condition (fT2l) is rearranged as (with k p = 1 and 

p « 1) 

2D -I v s R c R c {mi-m 2 ) ..,, 

m a > — a" (— ) = 2 (13) 

where m\ and rri2 are the inductor current slopes in the stages S\ and 52 respectively. Note that ([T3T l shows 
the required ramp slope m a to stabilize the subharmonic oscillation when the output capacitor is very large. The 
critical condition (TT3T > is reasonable, because, with very large output capacitor, the output voltage ripple is mostly 
contributed by the inductor current (multiplied by R c ) and the straight-line inductor current analysis (as in the 
current-mode control) is adequate for stability analysis. When the output capacitor is small, the critical condition 
(fT2l is more accurate, with an additional (second) term related to the output capacitor ripple. 
Without the ramp (m a = 0), rearranging (fl2l . subharmonic oscillation is avoided if 

T 1 

< — (14) 
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or, equivalently, 

R C C 1 D 2 1 

-|- > - + and D < - (15) 

T 2 1-2D 2 v ' 

also agreed with Eq. 5]. 
Remarks: 

(a) The conditions (fT2l is applicable to Z?o?/i PVMC and CF-PVR, and ( fT4l) or (fT51) is a special case (m a = 0) 
of CE2]). 

(b) L> < 1/2 is explicitly required in (fT5l . whereas I? < 1/2 is implicitly required in (fl~4T >. 

C PVMC wfrft fi c = 0: Case C 6 
Fro i? c = 0, the zero uj z is at infinity and Cj becomes Cq. Either from Table U or (fTTT i. the critical condition is 

^g ( 2^-2D + l) = l (16) 

agreed with Hi which claims that the ripple amplitude such as Av (or Ay) is an index to predict the subharmonic 
oscillation. In the buck converter, it happens that Av = (D 2 — D)T 2 v s /8LC HI, and the critical condition (f76l ) 
can be expressed as a condition in terms of Av . However, the condition (fT6l ) is valid only for R c = 0. 

D. PVMC with R c = and Small R: Case C 5 

As the load resistance R decreases, a pole at uj p = 1/ RC becomes significant. At high frequency, from (H), 

which is of case C5. The critical condition expressed in terms of v s is 

LV m ui s 



" s Rk p (a (D)-a(D,p)) (18) 

One sees that as R decreases, the critical value of v s increases proportionally according to l/R(a>o(D) — a(D,p)). 

E. PVMC with Simple RL Circuit: Case C3 

Most reported subharmonic oscillations occur in switched RLC circuits. However, subharmonic oscillation may 
occur even in a simple switched RL circuit of first order. 

Example 1. (Subharmonic oscillation occurs even with phase margin of 90.6° and infinite gain margin in the 
average model. ) Consider a simple RL circuit shown in Fig. [9] It is actually equivalent to peak current (voltage) 
control. From Fig. |6] for small p and D > 1/2, subharmonic oscillation may occur. 

Let the converter parameters be v s = 10 V, v r = 7.5, Vi = 0, Vh = 1, f s = 1 MHz, L = 1 /xH, and R = 1 Q. 

One has p = lo p /lo s = R/Llo s and 

T ^ = v 7TTJ ,m (19) 

V m (l + Ls/R) 

For k p = 9, the Bode plot (Fig. [T0l) shows a phase margin of 90.6° and an infinite gain margin. However, 
subharmonic oscillation still occurs as shown next. 

Based on the exact sampled-data model |[T2l . the pole as a function of k p is shown in Fig. \TT\ which shows the 
occurrence of subharmonic oscillation at k p = 8.63. Note that ( fT9l ) is of case C3 or C\. From Table U the critical 
condition is 

v s k p pa(D, p) = V m (20) 

Solving (l20l and the steady-state condition k p (v r — v ) = h(d), the critical gain is k p = 8.52, close to the exact 
sampled-data analysis. 

First, let k p = 8. The T-periodic orbit is stable as shown in Fig. [12] 

Next, let k p = 9. Subharmonic oscillation occurs with a 2T-periodic orbit as shown in Fig. [13] □ 
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Figure 9. Simple PVMC RL circuit. 
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Figure 10. Bode plot of the RL circuit shows a phase margin of 90.6° and an infinite gain margin. 



VIII. Average Current Mode Control (ACMC) with a Type-II Compensator: Case C 5 

The operation of ACMC is as follows |[T6l : The inductor current ii is sensed by a resistor R s and compared with 

a voltage reference v r from the voltage loop. The difference is amplified by a current-loop compensator, generally 

a type-II compensator, 

R S K C (1 + s/z c ) 



G c {s) 



(21) 



s(l + s/ujp) 
which has a small zero z c <C u) s , an integrator pole 0, a large pole uj p , and a gain K c . At high frequency, from (0, 

v s G c (s)G p (s) v s R s K c 1 



T{s) 



which is of case Cr, The critical condition is 
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(22) 



(23) 



It can be also expressed in terms of the ramp slope. From (I23t . the minimum ramp slope to avoid the subharmonic 
oscillation is 
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(24) 



A plot of ao(D) — a(D,p) (also a scaled L-plot C{D,p)) is shown in Fig. [T4j The corresponding contour plot 
is shown in Fig. [15] 
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Figure 11. Sampled-data pole as a function of k p . Subharmonic oscillation occurs at k v = 8.64. 





Figure 12. Stable T-periodic orbit, k p 



The contour plot is helpful to determine the maximum allowable loop gain without subharmonic oscillation. For 
example, suppose the ACMC converter is designed to operate at D = 0.7 and p = 0.4. From the contour plot, 
the L-plot value (cto(D) — a(D,p)) is around 1.5 (also the "elevation" or "contour level" in the contour plot). 
Let K = v s R s K c /V m z c Luj s . Then, from (|23l ), K < 2/3 is required to avoid the subharmonic oscillation. Since 
the largest value of ao(D) — a(D,p) is n, the subharmonic oscillation is completely avoided for any D or p if 
K < 1/vr. 

The contour plot also defines the critical boundary. For example, if K = 1, the contour plot with the elevation 
at 1 shown in Fig. [15] defines the critical boundary which separate stable region from the unstable region in the 
(D,p) space. Only those regions in the (D,p) space with the elevation less than 1 are the stable operating regions. 

The contour plot is also helpful to see the "weak spots" which are susceptible to subharmonic oscillation in 
the parameter space. Those regions with high elevations in the contour plot are the weak spots. In Fig. [15j one 
sees that for a large or small D, the elevation is high, and these regions are the weak spots. For mid-range of D, 
there exists a window of p such that the elevation is high. For example, in Fig. [151 draw a line at D = 0.6. As p 
increases, for p G (0.2,0.5), the plot rises above 1 and then falls below 1. In this window of p, the converter is 
prone to subharmonic oscillation. 
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Figure 13. Subharmonic oscillation with 2T-periodic orbit, k p — 9. 
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Figure 14. Plot of ao(D) — a(D,p) for case Cs. 



The instability window of p can be estimated. Since 

a (D) -a(D,p) « 

J a.\{D)p (for a small p) 

\ a {D) - 4vre~ 27rp + 2vre~ 27rDp (for a large p) 

Based on d23l ) and d25l) , the instability window of p is 

1 
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1 2 £>-l + 2e-^-^ . 
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(25) 



(26) 



l i^ai(Z))'2 

From (|23T ) and Fig. [15] for p being inside this window, with K ss 1 and 0.1 < D < 0.7, the subharmonic oscillation 
would occur, as shown in the next example. 



Example 2. {Instability window of p.) Consider a buck converter under ACMC 11171 p. 114]. The power stage 
parameters are v s = 14 V, v = 5 V, v r = 0.5, V L = 0, V h = 1, f s = 50 kHz, L = 46.1 /^H, C = 380 pF with ESR 
i? c = 0.02 Q, and i? = 1 Q,. The inductor current sensing resistance R s is 0.1 Q. The compensator has a zero at 
z c = 5652.9 rad/s, poles at and uj p , and a gain K c = 75506. 

Time-domain simulation. The compensator pole uj p is varied from 0.14w s to 0.81cj s . An unstable window of 
uj p between 0.18w s and 0.49w s was found and reported in JSJ. When w p is inside the window, the subharmonic 
oscillation occurs. 
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Figure 15. Contour plot of ao(-D) — a(D,p) for case C5. 
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Figure 16. Stable T-periodic solution, oj p = 0.15o> s 



For example, let co p = 0.15u; s . The converter is stable (Fig. HBV Let co p = 0A9co s . The converter is unstable with 
subharmonic oscillation (Fig. [17). Next, let lo p = 0.81w s . The converter is stable again (Fig. li~8V 

Independent sampled-data analysis. The sampled-data pole trajectories for uj p /uj s G (0.1,0.8) are shown in 
Fig. [19] There are four poles. Two poles are fixed around 0.88, and 0.95 (« e^c). A pole leaves the unit circle 
through -1 when ui p = 0.18w s , and enters the unit circle when uj p = 0A9to s . This explains exactly the instability 
window of uj p . 

Accurate prediction by the L-plot. The L-plot from (|23l i, as a function of p for D = v /v s = 0.357 is shown 
in Fig. [20] It shows the instability window of p € [0.18, 0.46]. The small error is due to the approximation of T(s) 
in (l22l ). The L-plot is also a scaled cross-section at D=0.357 in Fig. [T4j The instability window of p predicted by 
(l26l ) is [0.15,0.58], which is also a close estimation. 

In ACMC, the subharmonic oscillation is unrelated to the ripple amplitude (H. As uo p increases, the amplitude 
of y°(t) increases monotonously (because G c (s) is a low-pass filter). For uj p = 0A9u s , the converter is unstable 
with a small amplitude of y°(t), whereas for uj p = 0.81u; s , the converter is stable with a larger amplitude of y°(t). 
The ripple amplitude index to predict the subharmonic oscillation hypothesized in (6l does not apply in this case. 

In this example, the existence of instability window of to p is verified by time-domain simulation, pole trajectories, 
and agreed closely with the derived critical condition. □ 
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Figure 17. Subharmonic oscillation, uj p — 0.49o; s 




Figure 18. Stable T-periodic solution, uj p — 0.81oj s . 



IX. VMC with a Type-III Compensator: Case C 5 

A typical guideline ifTol, p. 412] popular in industry to set the parameters of the type-III three-pole-two-zero 
compensator is as follows. Set one pole at as an integrator, the second pole at 1/ ' R C C (to cancel the power stage 
zero), and the third pole at to p = oj s /2. Set the gain K c to adjust the phase margin and the crossover frequency. 
Set the two zeros at K z /y/LC and 1/yLC (to cancel the power stage poles), where k z is a zero scale factor to 
have additional flexibility to adjust the phase margin and the crossover frequency. The zero scale factor k z used 
in industry typically varies between 0.1 and 1.2. As will be shown later, a smaller value of k z may lead to the 
subharmonic oscillation. Taking into account the above guidelines, the compensator has a transfer function 

K C {1 + VICs/k z )(1 + y/LCs) 



G c {s) 



s(l + s/u p )(l + R c Cs) 



From CO, the loop gain is 



T{s) 



v s G c (s)G p (s) 
v s K c p 



1 



V m K z S(l + S/U)p) 

which is of case C5. The critical condition expressed in terms of v s is 

V m K z U} s 



(at high frequency) 



K c p(a (D)-a(D,p)) 
As discussed above, a smaller k z leads to a smaller critical source voltage v s 



(27) 

(28) 
(29) 

(30) 
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Figure 19. Sampled-data pole trajectories for ui p /u) s £ (0.1,0.8). 
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Figure 20. L-plot l !23t as a function of p. Its intersection with a horizontal line at 1 is the instability window of uj v . For C(p) > 1, 
subharmonic oscillation occurs. 



Example 3. (With phase margin of 38.9° based on the average model, the subharmonic oscillation still occurs.) 
Consider a buck converter with the type-Ill compensator (127T ). Exactly the same parameters as in the technical 
document 0J] are used: f s = 1/T = 300 kHz, L = 900 nH, C = 990 /jF, R = 0.4 £1, R c = 5 mtt, v r = 3.3 V, 
and V m = 1.5 V. For the compensator, K c = 7.78 x 10 4 , the zeros are l/2V~LC = 1.675 x 10 4 , 1/y/LC = 3.35 x 10 4 , 
and the poles are lo p = u s /2 = 9.425 x 10 5 , and 1/R C C = 2.02 x 10 5 . 

Time-domain simulation (Fig. |2~TI ) shows that the subharmonic oscillation occurs when v s = 16 V (D rs 0.206). 
This is also confirmed by the exact sampled-data analysis with a sampled-data pole at -1 when the subharmonic 
oscillation occurs. Based on the average model (with T(jui) as in d28l)), the loop gain T(juj) shown in Fig. [22] 
has a phase margin of 38.9° for v s = 16. The frequency response also shows an infinite gain margin because the 
phase never reaches -180°, which means that no matter how much v s increases (to increase the loop gain), the 
converter is expected to be stable based on the average model. However, the subharmonic oscillation still occurs 
when v s = 16. 

With D = 0.2, the prediction of the critical v s by (l30l ) is 17.1, close to the simulation result and the sampled-data 
analysis. As discussed above, the small error is due to the approximation of T(s) in ( [29l . □ 

Since the VMC buck converter with a type-Ill compensator is of case C5, same as the ACMC buck converter. 
The VMC buck converter with a type-Ill compensator also has an instability window of ui p as shown in the next 
example. 
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Figure 21. Signal waveforms showing the subharmonic oscillation. 



Gm = Inf dB (al Inf rad/sec| , Pm = 38.9 deg (at 1 .08e+006 rad/sec) 




Frequency (rad/sec| 

Figure 22. Loop gain frequency response shows a phase margin of 38.9° and an infinite gain margin, but the subharmonic oscillation still 
occurs. 



Example 4. (Unstable window of uj p , unrelated to the ripple amplitude of y°(t).) Consider again Example [3] with 
v s = 16. Vary u p from 0.1w s to 0.6o; s . An instability window of uj p G (0.23, 0.5)w s is found. Similar to Example[2] 
the value of u p adjusts the ripple amplitude of y°(t). A larger uj p leads to a larger ripple of y°(t). In 151. it is 
hypothesized that the ripple amplitude of y°(t) is related to subharmonic oscillation. The following simulation 
shows that the ripple amplitude of y°(t) is unrelated to subharmonic oscillation even for VMC. 

Time-domain simulation. For uj p = 0.2uj s , the ripple amplitude of y°(t) is small, and the converter is stable 
(Fig. [23]). For uj p = 0.24w s , the ripple amplitude of y°(t) is larger, and the converter is unstable with subharmonic 
oscillation (Fig. l24b . For to p = 0.6w s , the ripple amplitude of y°(t) (not shown) is even larger, but the converter 



is stable (Fig. I25l ). Comparing Figs. 12311251 the ripple amplitude of y°(t) is unrelated to subharmonic oscillation. 
This also shows a counter-example for the hypothesis proposed in [6] that the ripple amplitude of y°(t) is related 
to the subharmonic oscillation. 

Independent sampled-data analysis. The instability window of uj p is also confirmed by the sampled-data pole 
trajectories. The sampled-data pole trajectories for 0.1w s < lo p < 0.6uj s are shown in Fig. |26j There are five poles. 
Three poles are fixed around 0.9485, 0.8853, and 0.51. A pole leaves the unit circle through -1 when u p = 0.23o; s , 
and enters the unit circle when uj p = 0.5u s . This explains exactly the instability window of ui p . 

Accurate prediction by the L-plot. The L-plot (where T(s) is based on (|29l ) is shown in Fig. |27]which indicates 
an instability window of uj p G [0.23, 0A7]u s . As discussed above, the small error is due to the approximation of 
T(s) in (|29l . The instability window of p predicted by (l26l i is [0.17,0.58], which is also a close estimation. 

In this example, the existence of instability window of u p is also verified by time-domain simulation, pole 
trajectories, and agreed closely with the derived critical condition. □ 
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Figure 23. The converter is stable, uj p = 0.2uj s 
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Figure 24. Subharmonic oscillation, u p — 0.24o; s 



X. Conclusion 

A general critical condition d4) of subharmonic oscillation in terms of the loop gain is derived. It is applicable to 
a general nonlinear switching system represented in Fig. [5] Therefore, it is also applicable to other similar nonlinear 
control systems. Many closed-form critical conditions for various control schemes in terms of converter parameters 
are also derived. They are summarized in Table |II] The effects of different parameters (such as v s , R, R c , and the 
ramp slope m a ) on the subharmonic oscillation can be clearly seen. 

The questions asked in the Introduction are answered: 

1) Those previously known critical conditions, such as (J9j for CMC, (IT2T > for V 2 control, and (fT6l ) for VMC 
with zero ESR, become special cases in the generalized framework. 

2) The instability window of pole can be explained by case C5. A system of the case C5, such as ACMC or 
VMC with a type-Ill compensator, will have an instability window of pole. 

3) The hypothesis [6] that the signal ripple amplitude can predict the occurrence of subharmonic oscillation is 
applicable for a converter with zero ESR (case Cq). It is not applicable to a converter with a large pole (case 
C5). For example, in ACMC or VMC with a type-Ill compensator, the instability window of pole indicates 
that the signal ripple amplitude is not an accurate index to predict subharmonic oscillation. Also note that the 
ripple amplitude, either Av or Ail for R c = 0, has a term of D(l — D) which is symmetric with respect 
to D = 1/2. For a ripple index to be valid, the critical condition also needs to be symmetric with respect 
to D = 1/2. Among the cases in Table U only the critical condition for Cq is symmetric with respect to 
D = 1/2. Therefore, the ripple index has very limited applications. 

4) A typical critical condition is a weighted combination of three terms: ao(D) = tt(2D — 1), a\{D) = 
7r 2 (2D 2 — 2D + 1), and a correction term c(D,p). For example, both V 2 control and PVMC have the same 
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Figure 25. The converter is stable with a larger ripple of y°(t), uj p = 0.6uj s 
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Figure 26. Sampled-data pole trajectories for 0.1uu 3 < uj p < 0.6lo 3 



fonn of critical condition. Also, both PVMC and CMC have a term ao(D) = ir(2D — 1). For a compensator 
with a pole smaller than one tenth of the switching frequency, the correction term can be ignored. 
5) Given an arbitrary control scheme, a systematic procedure is proposed to derive the critical condition for 
that control scheme. First, approximate the loop gain T(s) at high frequency (higher than ui s /2). Then, from 
Table H one can readily obtain the critical condition in terms of converter parameters. 
Given the closed-form condition in terms of converter parameters, one knows the quantitative effect of each 
parameter on the subharmonic oscillation. One can make the L-plot (such as Figs. [20] and |2"71 ) as a function of 
a parameter of interest, its intersection with a horizontal line at 1 determines the stable operating range of that 
parameter. Also, given parameters pi and p2 for example, based on the closed-form condition, one can determine the 
instability boundary in the parameter space (pi, P2) as shown in Fig. [Q or from the contour plot (such as Fig. H3V 
Based on these plots, one knows how to choose proper parameter values to avoid the subharmonic oscillation. 

This paper focuses on the converter operated at a fixed switching frequency. Similar analysis can be applied to 
the converter with variable frequency control, such as constant on-time control. The F-transforms in Table I still 
apply, but with a different a(D,p). The results are reported separately. 



References 

[1] R. W. Erickson and D. Maksimovic, Fundamentals of Power Electronics, 2nd ed. Berlin, Germany: Springer, 2001. 

[2] D. C. Hamill, J. H. B. Deane, and J. Jefferies, "Modeling of chaotic DC-DC converters by iterated nonlinear mappings," IEEE Trans. 

Power Electron., vol. 7, no. 1, pp. 25-36, 1992. 
[3] R. B. Ridley, "A new, continuous-time model for current-mode control," IEEE Trans. Power Electron., vol. 6, no. 2, pp. 271-280, 1991. 
[4] F. D. Tan and R. D. Middlebrook, "A unified model for current-programmed converters," IEEE Trans. Power Electron., vol. 10, no. 4, 

pp. 397^08, 1995. 



18 




0.4 0.6 



Figure 27. L-plot as a function of p. Its intersection with a horizontal line at 1 is the instability window of uj p . For C(p) > 1, subharmonic 
oscillation occurs. 

Table II 
Critical conditions for different control schemes. 



C 2 CMC, %(D 



-n 



v s k pP T* r2R a C , 



C 7 PVMC, R c / 0, "tv^lc [^^(2D - 1) + (2D 2 - 2D + 1) 
C 7 V 2 Control, same as PVMC above 

2 2D + 1) = 1 



= 1 



C 6 PVMC with R c = 0, lt P J n (2D 



4V m LC* 

C s PVMC with R c = and Small 7?, v s = 

C 3 PVMC RL circuit, « s A;ppa(D,p) = V, 

Ch ACMC with type-II compensator, v 3 — 

Ch VMC with type-Ill compensator, v 3 = 





LV m ^ 








Rk p 
-n. 


("11(D)- 
V m z c Li 


a 


D 


P)) 


R S K C 


(«o(D)- 
Vrn k-z^s 


-CL 


D 


.;.)) 



K c p(a (D)-a(D,p)) 



[5] C.-C. Fang, "Sampled data poles, zeros, and modeling for current mode control," International Journal of Circuit Theory and 

Applications, 2011, accepted and published online, DOI: 10. 1002/cta.790. 
[6] A. El Aroudi, E. Rodriguez, R. Leyva, and E. Alarcon, "A design-oriented combined approach for bifurcation prediction in switched- 

mode power converters," IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 57, no. 3, pp. 218-222, Mar. 2010. 
[7] R. Redl and J. Sun, "Ripple-based control of switching regulators - an overview," IEEE Trans. Power Electron., vol. 24, no. 12, pp. 

2669-2680, Dec. 2009. 
[8] C.-C. Fang, "Modeling and instability of average current control," in EPE International Power Electronics And Motion Control 

Conference, 2002, paper SSIN-03. 
[9] J. Li and F. C. Lee, "Modeling of V 2 current-mode control," IEEE Transactions on Circuits and Systems-l: Regular Papers, vol. 57, 

no. 9, pp. 2552-2563, 2010. 
[10] R. Genesio and A. Tesi, "Harmonic balance methods for the analysis of chaotic dynamics in nonlinear systems," Automatica, vol. 28, 

no. 3, pp. 531-548, 1992. 
[11] A. Tesi, E. H. Abed, R. Genesio, and H. O. Wang, "Harmonic balance analysis of period-doubling bifurcations with implications for 

control of nonlinear dynamics," Automatica, vol. 32, no. 9, pp. 1255-1271, 1996. 
[12] C.-C. Fang, "Sampled-data analysis and control of DC-DC switching converters," Ph.D. dissertation, Dept. of Elect. Eng., Univ. of 

Maryland, College Park, 1997, available: http://www.lib.umd.edu/drum/ also published by UMI Dissertation Publishing in 1997. 
[13] C.-C. Fang and E. H. Abed, "Harmonic balance analysis and control of period doubling bifurcation in buck converters," in Proc. IEEE 

ISCAS, vol. 3, May 2001, pp. 209-212. 
[14] , "Analysis and control of period-doubling bifurcation in buck converters using harmonic balance," Latin American Applied 

Research: An International Journal, pp. 149-156, 2001, Special theme issue: Bifurcation Control: Methodologies and Applications, In 

Honor of the 65th Birthday of Professor Leon O. Chua. 
[15] M. K. Kazimierczuk, "Transfer function of current modulator in PWM converters with current-mode control," IEEE Transactions on 

Circuits and Systems-l: Fundamental Theory and Applications, vol. 47, no. 9, pp. 1407-1412, 2000. 
[16] C. P. Basso, Switch-Mode Power Supplies. McGraw-Hill, 2008. 
[17] W. Tang, F. C. Lee, and R. B. Ridley, "Small-signal modeling of average current-mode control," IEEE Trans. Power Electron., vol. 8, 

no. 2, pp. 112-119, 1993. 



19 

[18] D. Mattingly, "Designing stable compensation networks for single phase voltage mode buck regulators," Intersil Americas Inc., Tech. 
Rep., 2003, available: www.intersil.com/data/tb/tb417.pdf. 



